Peregrination of the selectivity filter delineates the pore of the human 
voltage-gated proton channel hH v 1 



Deri Morgan, 1 Boris Musset, 1 Kethika Kulleperuma, 2,3 Susan M.E. Smith, 4 Sindhu Rajan, 5 
Vladimir V. Cherny, 1 Regis Pomes, 2,3 and Thomas E. DeCoursey 1 

'Department of Molecular Biophysics and Physiology, Rush University, Chicago, IL 60612 
2 Molecular Structure and Function, Hospital for Sick Children, Toronto, Ontario M5G 1X8, Canada 
3 Department of Biochemistry, University of Toronto, Toronto, Ontario M5S 1A1, Canada 
"Department of Biology and Physics, Kennesaw State University, Kennesaw, GA 30144 
5 Department of Medicine, University of Chicago, Chicago, IL 60637 

Extraordinary selectivity is crucial to all proton-conducting molecules, including the human voltage-gated proton 
channel (hH v l), because the proton concentration is >10 6 times lower than that of other cations. Here we use 
"selectivity filter scanning" to elucidate the molecular requirements for proton-specific conduction in hH v l. Asp 112 , 
in the middle of the SI transmembrane helix, is an essential part of the selectivity filter in wild-type (WT) channels. 
After neutralizing Asp 112 by mutating it to Ala (D112A), we introduced Asp at each position along SI from 108 to 
118, searching for "second site suppressor" activity. Surprisingly, most mutants lacked even the anion conduction 
exhibited by D112A. Proton-specific conduction was restored only with Asp or Glu at position 116. The D112V/ 
VI 16D channel strikingly resembled WT in selectivity, kinetics, and ApH-dependent gating. The S4 segment of this 
mutant has similar accessibility to WT in open channels, because R211H/D112V/V116D was inhibited by inter- 
nally applied Zn 2+ . Asp at position 109 allowed anion permeation in combination with D112A but did not rescue 
function in the nonconducting D112V mutant, indicating that selectivity is established externally to the constric- 
tion at F150. The three positions that permitted conduction all line the pore in our homology model, clearly de- 
lineating the conduction pathway. Evidently, a carboxyl group must face the pore directly to enable conduction. 
Molecular dynamics simulations indicate reorganization of hydrogen bond networks in the external vestibule in 
D112V/V116D. At both positions where it produces proton selectivity, Asp frequently engages in salt linkage with 
one or more Arg residues from S4. Surprisingly, mean hydration profiles were similar in proton-selective, anion- 
permeable, and nonconducting constructs. That the selectivity filter functions in a new location helps to define 
local environmental features required to produce proton-selective conduction. 



INTRODUCTION 

Voltage-gated proton channels (H v ls) enable phago- 
cytes to kill pathogens (Henderson et al., 1988; DeCoursey 
et al., 2003; Morgan et al., 2009; DeCoursey, 2010; 
Demaurex, 2012), basophils to secrete histamine (Musset 
et al., 2008b), and airway epithelia to control surface pH 
(Fischer, 2012), as well as enable sperm motility (Musset 
et al., 2012) and capacitation (Lishko et al., 2010), and 
B lymphocyte signaling (Capasso et al., 2010), and may 
exacerbate breast cancer metastasis (Wang et al., 2012) 
and ischemic brain damage (Wu et al., 2012). All of 
these functions are predicated on the proton specificity 
of H v l. The low concentration of H + in biological fluids 
means that extraordinary selectivity is necessary even to 
ensure that H + is the main conducted species. In fact, 
proton selectivity in H v l appears to be perfect (Musset 
et al., 2011; DeCoursey, 2013). 
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An acidic group in the middle of the SI transmem- 
brane segment is critical to the proton specificity of H v l 
and is provided by Asp 112 in human H v l (hH v l; Musset 
et al., 2011) and Asp 51 in H v l from a dinoflagellate, Kar- 
lodinium veneficum (Smith et al., 2011). Despite only 15% 
amino acid identity of the proteins, the conservative 
Asp— >Glu mutation preserved proton specificity, whereas 
Ser, Ala, or His substitution for Asp at this position re- 
sulted in anion permeability in both species, strongly 
suggesting that the selectivity mechanism is widely con- 
served evolutionarily. The presence of an Asp facing 
the pore is not sufficient, however, because Asp 185 can 
be neutralized without compromising proton selectiv- 
ity, and does not preserve selectivity when Asp 112 is neu- 
tralized (Musset et al., 2011). Other molecular elements 
that may be required are not known. Our homology 
model indicates that the second of three Arg residues in 
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the S4 segment, Arg , forms a salt bridge with Asp , 
and that the resulting charge compensation is impor- 
tant for proton selectivity (Kulleperuma etal., 2013). To 
refine further the molecular requirements of the selec- 
tivity filter, we explore here the extent to which the criti- 
cal Asp can be moved along the SI segment. We find 
that an excellent proton channel is produced when Asp 
is shifted from position 1 12 to position 116. The mutant 
channel is proton specific, exhibits ApH-dependent gat- 
ing characteristic of all H v ls, and surprisingly, neutral 
amino acid substituents at this location produce anion 
permeability. Molecular dynamics (MD) simulations in- 
dicate that Asp 116 forms a salt bridge with Arg 205 and/ or 
Arg 208 , but the latter also pairs with Asp 183 , reflecting 
reorganization of charge clusters in the mutant com- 
pared with WT channels. The results underline the im- 
portance of intramolecular charge compensation for 
proton selectivity. 

MATERIALS AND METHODS 

Gene expression 

Site-directed mutants were created using the Stratagene QuikChange 
(Agilent Technologies) procedure according to the manufactur- 
er's instructions. Transfection was done as described previously 
(Kulleperuma et al., 2013). Both HEK-293 cells and COS-7 cells 
were used as expression systems, the latter more frequently. In 
a previous study, we systematically compared the properties of 
liH v l when expressed in these two cell lines and found no differ- 
ence (Musset et al., 2008a). Although currents that decayed at 
large positive voltages (presumed to be volume-regulated anion 
currents) were sometimes seen at the start of experiments, these 
disappeared over time. Occasional cells displayed a few unidenti- 
fied single-channel currents superimposed on the macroscopic 
currents. The unitary conductance of hH v l is just 140 fS at pH; 
5.5 (Cherny et al., 2003); thus, visible unitary currents were con- 
sidered extraneous and were ignored. No other voltage- or time- 
dependent conductances were observed under the conditions of 
this study. Thus, depolarization-activated time-dependent cur- 
rents, which in many cells were orders of magnitude larger than 
any background currents, were assumed to reflect the transfected 
construct. Both cell lines sometimes exhibit small native H v l cur- 
rents, which could be distinguished from transfected mutant 
channels by their high Zn 2+ sensitivity (see Results). 

Electrophysiology 

GFP-tagged proton channels were identified using inverted mi- 
croscopes (Nikon) with fluorescence capability. Conventional 
patch-clamp techniques were used at 21°C or at room tempera- 
ture (20-25°C) (Kulleperuma et al., 2013) . Bath and pipette solu- 
tions contained 60-100 mM of buffer, 1-2 mM CaCl 2 or MgCl 2 
(intracellular [IC] solutions were Ca 2+ free), 1-2 mM EGTA, and 
TMAMeS0 3 to adjust the osmolality to ~300 mOsm, titrated with 
TMAOH. Buffers used were Mes at pH 5.5-6.0, HomoPIPES at 
pH 4.5, and PIPES at pH 7.0. For Zn 2+ measurements, EGTA was 
omitted. Currents are shown without leak correction. Reversal po- 
tentials were corrected for measured liquid junction potentials. 
Unless stated otherwise, cells were held at a holding potential, 
Kola, before pulses and returned to Vh 0 i d after families of pulses. 

Reversal potentials were determined by two methods, depend- 
ing on the relative positions of Threshold and V ler For most con- 
structs, Vth r e S hoid was positive to V,„, and the latter was determined 



by examining tail currents (e.g., Fig. 2 B). Because hH v l currents 
were the only time-dependent conductance present, estimates of 
the amplitude and direction of current decay during deactivation 
were used to establish V rev (Morgan and DeCoursey, 2007) . By this 
procedure, time-independent leak or other extraneous conduc- 
tances do not affect V r „. For mutants in which Threshold was nega- 
tive to 14ev (for example, D112A/V116D in Fig. 2 C), it was possible 
to observe the reversal of the direction of currents activated dur- 
ing pulse families. Tail currents were not observed in nontrans- 
fected cells; for example, Fig. 1 B illustrates the absence of tail 
currents in a cell with the nonconducting D112V mutant. 

MD simulations 

MD simulations of the WT protein, single-point mutants D112V 
and D112S, and double-point mutants D112V/V116D ("VD") and 
D112V/V112S ("VS") were performed in a hydrated lipid bilayer 
based on the homology model constructed and validated in a re- 
cent study (Kulleperuma et al., 2013). 

12 conformations of the WT protein with pore-associated water 
were used as initial structures. These snapshots correspond to the 
endpoint of 12 different 200-ns-long unrestrained simulations 
in a membrane-mimetic octane slab (Kulleperuma et al., 2013). 
Each conformation represents one of the three configurations of 
the D112-R208 salt bridge: bidentate and monodentate confor- 
mations (involving two or one hydrogen bond, respectively) , and 
open, in order of increasing separation. A preequilibrated con- 
figuration of a l-palmitoyl-2-oleoylphosphatidylcholine (POPC) 
bilayer was obtained from a previous study (Kulleperuma et al., 
2013). The OPLS-AA protein force field (Jorgensen et al., 1996) 
was mixed with the Berger lipid parameters (Berger et al., 1997) 
by applying the half-e double-pairlist method (Chakrabarti et al., 
2010). The TIP3P water model was used (Jorgensen et al., 1983). 
InflateGRO (Kandt et al., 2007) was used to embed the protein in 
the bilayer. The system was then hydrated, and 54 Na + and 56 CP 
ions were added to neutralize the charge of the system and yield 
an approximate ionic concentration of 500 mM. The resulting 
simulator cell consisted of 126 POPC and ~5,900 water molecules 
in a box of ^6.5 x 6.5 x 8 nm 3 . The MD parameters used for this 
study are described elsewhere (Kulleperuma et al., 2013). 

Each of the 12 WT systems was first energy-minimized using 
50,000 steps of steepest descent, followed by an equilibration phase 
of 50 ns with position restraints on protein backbone and pore- 
associated water oxygen atoms. The production run consisted 
of 200-ns-long unrestrained simulations for each system. Snap- 
shots of protein and pore-associated water molecules were se- 
lected from each of the 12 equilibrated WT protein systems at 
t= 100 ns to produce single- and double-point mutants. Mutations 
were introduced using an in-house script, followed by 1,000 steps 
of energy minimization. Asp 112 and Val 116 side chains were modi- 
fied to Val, Ser, or Glu by either overwriting the WT heavy atoms 
or deleting some or both. Side-chain and backbone dihedral 
angles were checked after energy minimization. Extra ions were 
added to the solution as required to neutralize the system after 
mutations were introduced. Another 50,000 steps of energy mini- 
mization were performed before an equilibration phase consist- 
ing of an additional 25 ns with position restraints, as described 
above. 12 time trajectories of 200 ns differing in the initial confor- 
mation of the protein were generated for each mutant. The total 
production time was 12 us. For the final analysis, eight replicas of 
WT and VD and seven replicas of VS, D 1 1 2V, and D 1 1 2S were se- 
lected after discarding replicas that displayed significant changes 
in secondary structure. 

Snapshots saved every 20 ps during the last 100 ns of each 
selected production run were analyzed for each system. Molecu- 
lar graphics were generated by VMD 1.8.7 (Humphrey et al., 
1996), and all trajectories were analyzed using Gromacs tools and 
in-house codes. 
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Online supplemental material 

Figs. SI and S2 illustrate the structural plasticity of the channel 
with respect to translation of the four transmembrane helices 
during extended MD simulations. Table SI provides V lm data for 
D112V/V116D and D112V/V1 16D/R211H mutants in the pres- 
ence of several cations or Cl~ that show both constructs to be proton 
selective. The online supplemental material is available at http:// 
www.jgp.org/cgi/content/full/jgp.201311045/DCl. 

RESULTS 

We generated a series of mutants in which the Asp resi- 
due critical to proton specificity was effectively shifted 
up and down the SI helix to each position from 108 to 
118. In the initial series of experiments, we replaced Asp 112 
with Ala (D112A) to produce an anion-permeable chan- 
nel (Musset et al., 2011), and then introduced a second 
mutation with the goal of restoring proton selectivity 
("second site suppression") . Additional studies were done 
in the D112V background, which in a sense is more rig- 
orous, because this single mutant does not conduct at 
all (Musset et al., 201 1) . All mutations were introduced 
into a Zn 2+ -resistant background (H140A/H193A), so 
that spurious small native H v l currents often present in 
COS-7 or HEK-293 cells (Musset et al., 2011) could be 
identified by their sensitivity to 10 pM Zn 2+ . Distinct 
Zn 2+ -insensitive currents were observed only in mu- 
tants with Asp at positions 116 or 109 (D112A/V116D, 
D112A/V109D). 

The second site mutation V1 16D restores proton 
selectivity to nonconducting or anion-permeable mutants 
Fig. 1 shows that although the single-point mutation 
D112V eliminates current altogether (Musset et al., 2011) , 
introducing Asp at position 116 restores robust proton 
current to the double mutant, D112V/V116D. Similarly, 
the D112A single mutant is anion permeable (Musset 
et al., 201 1 ) , but introducing Asp at position 116 restored 
proton-specific current to the double mutant, D112A/ 
V116D. The proton selectivity of both double mutants, 
D112A/V116D and D112V/V116D, was confirmed by 
the proximity of their reversal potentials, V rev , to the 
Nernst potential for H + , E H (Fig. 2 A), over a wide pH 
range (pH D 4.5-7.5 and pH ; 5.5-7.0) . Fig. 2 B illustrates 
determination of V rev from tail currents in a cell with 
D112V/V116D channels at pH G 5.5 and 7.0, with pH ; 
5.5. As indicated by the arrows, V rev shifted from 0 to 
-77 mV, near E H of -87 mV. The D112A/V116D mu- 
tant activated in a more negative voltage range so that 
inward currents were observed negative to £h in fami- 
lies of currents. Fig. 2 C illustrates currents from pulses 
bracketing V rev that reveal an ^60-mV shift between 
pH Q 5.0 and 6.0. In addition, substituting Na + , Li + , or K + 
for TMA + , or CT for CH3SO3 had no effect on ^ ev 
(Table SI). In summary, shifting Asp from position 112 
to 1 16 moves the proton selectivity filter outward by one 
turn of the helix. 



Asp supports current only when facing the pore 

Surprisingly, moving Asp to positions other than 109 or 
116 eliminated voltage-gated current altogether. All mu- 
tants were tagged with GFP, and transfected cells with 
green "halos" indicating membrane expression were se- 
lected under fluorescence for recording. Nonconduct- 
ing Asp mutants included (mutation, number of cells) : 
positions 108 (D112A/L108D, 10), 110 (D112A/V110D,4), 
111 (D112A/L111D, 4), 113 (D112A/A113D, 5), 114 
(D112A/L114D,4), 115 (D112A/L115D, 4), 117 (D112A/ 
L117D, 7), and 118 (D112A/A118D, 4). Positions 116, 
112, and 109 (Fig. 1, cartoon) all face the "pore" in the 
predicted open-state structure of hH v l (Kulleperuma 
et al., 2013). Evidently, Asp at a nonpore-facing location 
fails to support conductance, confirming that positions 
109, 112, and 116 line the pore of hH v l. 

As a control for the possibility that nonconducting 
mutants did not fold properly, we created a single mu- 
tant, Al 13D, in which the native Asp 112 was preserved. If 
Asp at a nonpore-facing location caused the protein to 
misfold, A113D should not conduct. Instead, A113D 
displayed small voltage-gated currents that were kineti- 
cally different from WT currents but were unequivocally 




Figure 1. The D112V mutation abolishes current, but V116D 
restores proton-specific current, an example of second-site sup- 
pression. Whole-cell currents at pH Q 7.0 and pH; 5.5 during pulses 
in 10-mV increments up to the indicated voltages for WT (A), 
D112V (B), or D112V/V116D (C), all expressed in COS-7 cells. 
Holding potential ( Vhoid) and pulse durations were —90 mV, 1 s 
(A), -40 mV, 3 s (B), and -60 mV, 2 s (C). In all cases, the volt- 
age was returned to V^oid after the pulse, which is why the tail cur- 
rents are inward for A and outward for C. Cartoons in all figures 
indicate SI and S4 helices, with color coding as follows: red, Asp 
or Glu; yellow, Val; blue, Arg; gray, other amino acids or nonpore- 
facing residues. 
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proton selective (Fig. 3). Thus, nonconducting mu- 
tants, including D112A/A113D, most likely were ex- 
pressed but nonfunctional. If, despite the appearance 
of green protein in the membrane, some mutants mis- 
folded, the fact remains that these proteins do not func- 
tion as channels. 

At position 109, Asp plays a permissive role 

Introducing Asp at position 109 into the nonconduct- 
ing D112V background (Musset et al., 2011) did not 
overcome the lack of conductance produced by the 
D 1 1 2V single mutation ( n = 8 cells) . However, the D 1 1 2A/ 
V109D mutant exhibited distinct currents at pH Q 5.5 
and pH ; 5.5 (Fig. 4 A), in contrast to the majority of 
mutants that did not conduct (D112A/L108D, D112A/ 
V110D, D112A/L111D, D112A/A113D, D112A/L114D, 
D112A/L115D, D112A/L117D, and D112A/A118D). 
Replacing CH 3 S03 with CP increased the outward 
current (Fig. 4 C) , reflecting CP influx, and pro- 
duced a large negative shift of V m (Fig. 4, B vs. D), con- 
firming CP permeability. The shift of V rev when CP 
replaced CH 3 S0 3 was -35.9 + 4.2 mV (SEM; n = 5), in 
the range reported for the Dl 12A single mutant, —29 mV 



(Musset et al., 201 1) . Evidently, introducing Asp at posi- 
tion 109 did not interfere with anion conduction seen 
in the D112A single mutant (Musset et al., 2011). Both 
the anion selectivity of D112A/V109D and the lack of 
conductance in D112V/V109D suggest that at position 
109, Asp cannot mediate proton selectivity and has no 
discernible effect on the selectivity that is established 
elsewhere. The two positions where Asp produced proton 
selectivity are in the external vestibule in our model, 
outside the highly conserved charge transfer center de- 
limiter Phe 150 (Tao et al., 2010). 

The engineered Asp 116 proton channel functionally 
resembles WT 

To what extent does moving the selectivity filter out- 
ward by one turn of the helix reinstate native hH v l 
properties? The gating and pH dependence of the 
D112A/V116D and D112V/V116D mutants are illus- 
trated in Fig. 5. The general appearance of the proton 
currents is unremarkable; the voltage dependence and 
gating kinetics are roughly similar to WT. In D112V/ 
V116D, the time constant of tail current decay, r tai i, was 
1.1 + 0.02 s (mean + SEM; n = 3) at -40 mV and pH 




Figure 2. Second site suppression in double mu- 
tants. (A) Shifting the crucial aspartate from posi- 
tion 112 to 116 (V116D) restores proton selectivity 
to both the nonconducting D112V (red symbols) 
and the anion-permeable D112A (blue symbols) 
single mutants. Measurements in the same cell are 
connected by lines. For both, V t „ measured over a 
wide range (for D112V/V116D pH 0 was 5.5, 6.0, 6.5, 
7.0, 7.5, or 8.0, and pHj was 5.5, 6.0, 6.5, or 7.0; for 
D112A/V116D pH 0 was 4.5, 5.5, 6.0, 6.5, or 7.0, and 
pH; was 5.5, 6.0, or 6.5) falls close to the Nernst po- 
tential for H + , £ H (dashed green line). (B) Determi- 
nation of V la . from tail currents for Dl 12V/V1 16D in 
a COS-7 cell with pH ; 5.5. Prepulses to +30 mV, pH 0 
5.5, or to —40 mV, pH„ 7.0, activated the conduc- 
tance, followed by repolarization to the indicated 
voltages in 10-mV increments, with the most posi- 
tive labeled. (C) Determination of V IC , from current 
families for D112A/V116D in a COS-7 cell with pH ; 
5.5. As indicated, V ho id was -60 mV (left) or -70 mV 
(right); the voltage was returned to Vh old after the 
pulses. Currents are shown during selected pulses 
that bracket V r „ in 10-mV increments. Inward cur- 
rent was activated negative to V T „, and the first out- 
ward current is labeled. 
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5.5//5.5, which is close to the WT value of 0.81 s mea- 
sured under the same conditions (Musset et al., 2011). 
The activation time constant, T act , was 2.3 ± 0.3 s (n = 3) 
at +40 mV in D112V/V116D, compared with the WT 
value of 1.1 s (Musset et al., 2011). Thus, marked 



Figure 3. Introducing Asp at position 113, pre- 
dicted by our model to be in a nonpore-facing 
location, results in membrane expression of a func- 
tioning proton-selective channel. The cartoon 
emphasizes that Asp 112 is still present. Families 
of currents are shown in a COS-7 cell in whole- 
cell configuration at pH c 5.5 (A) or 7.0 (B), with 
pH; 5.5, with pulses applied from as labeled 
to the indicated voltages in 10-mV increments. 
Cells were returned to V boM after pulses. (C) Pro- 
ton selectivity is shown by the proximity of V T „ 
to -Eh (dashed line). Insets show V ra , determina- 
tion (left) at pH 0 5.5 and pHj 5.5 by reversal of 
current during a family of pulses in 10-mV incre- 
ments (VhoM = —60 mV) and at pH„ 7.0 and pH ; 5.5 
(right) by tail currents. V r „ was measured at pH 0 
4.5, 5.0, 5.5, 6.5, 7.0, and 7.5, and at pH, 5.5 or 6.5. 

changes in the gating of the mutant channels were 
not observed. 

A characteristic property shared by all known H v l 
is tight regulation of the position of the gu~V relation- 
ship by the pH gradient, ApH (Cherny et al., 1995; 
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Figure 4. Moving aspartate from po- 
sition 112 to 109 results in anion cur- 
rents. Whole-cell currents in a COS-7 
cell expressing D112A/V109D chan- 
nels, all at pH Q 5.5 and pH ; 5.5, in sym- 
metrical TMA + CH 3 SCV (A and B) or 
with Cr in the bath (C and D). Pulses 
applied in 10-mV increments. V hold 
was -40 mV (A and B) or -60 mV 
(C and D) . Cells were returned to Vh 0 i d 
after pulses. V ICY determination from 
tail currents (B and D), with pulses in 
10-mV increments. 
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DeCoursey, 2013). Because the mechanism remains mys- 
terious, we examined this property. As evident in Fig. 5 
(D and E) , changes in pH Q produced shifts in both In-V 
and gH~V curves of roughly 40 mV/U, similar to the 
shifts observed in native and WT hH v ls (Ramsey et al., 
2006; Musset et al., 2008a). Plots of Threshold versus V rev 
for D112V/V116D (Fig. 5, D and E, red) and D112A/ 
V116D (blue) had slopes (0.64 and 0.63, respectively) 
similar to those observed previously for WT hH v ls 
(Musset et al., 2008a). The D112A/V116D mutant acti- 
vated at voltages 26 mV more negative than D112V/ 
V116D, suggesting relative stabilization of the open 
state by Ala versus Val at position 1 1 2. As a consequence 
of this shift, threshold in D112A/V116D was often nega- 
tive to £h so that inward currents were observed. Never- 
theless, ApH-dependent gating was preserved in both 
double mutants. 



Mutations at position 116 (in D112V channels) mimic 
mutations at 112 

That proton selectivity was restored to D112V and 
Dl 12A mutants by introducing Asp at position 116 sug- 
gests that the intramolecular interactions at both posi- 
tions that contribute to proton selectivity are similar. To 
explore the extent of equivalence of these positions, we 
compared effects of several other mutations at 1 16 with 
those at 112. Fig. 6 shows that robust proton-selective 
currents were observed when Glu replaced Asp at posi- 
tion 116 (in the D112V background), just as with re- 
placement of Asp by Glu at position 112 (Musset et al., 
2011). We found previously that neutral mutants of Asp 112 
were anion selective (Musset et al., 201 1) . Astoundingly, 
both Asn and Ser at position 116 produced anion-per- 
meable channels in the D112V background, which it- 
self does not conduct. Fig. 7 illustrates a current family 




Figure 5. Voltage and pH dependence of gating of hH v l mutants with the selectivity filter shifter from 112 to 116. Families of cur- 
rents in a COS-7 cell expressing D112A/V116D channels at pH ; 6.5 and pH Q 7.5 (A), 6.5 (B), or 5.5 (C) are shown in 10-mV increments 
as labeled, from 14ou — 60 mV (A) or —30 mV (B and C). Cells were returned to Vhoid after pulses. (D) Current-voltage relationships 
in this cell. Currents were fitted to a rising exponential function and extrapolated to infinite time. Note inward currents at pH„ 5.5. 
(E) Conductance-voltage relationships for the same currents. Limiting slope conductance for the most negative voltages provides a gat- 
ing charge estimate of 5.3-6.0 (F) Regulation of voltage gating by ApH. Ktoshoid is plotted against V rev measured in the same cell and 
solution. Lines show linear regression fits defined by threshold = 0.64 t^ ev + 0 mV (Dl 12V/V1 16D; red circles) or Kh resllol( i = 0.63 V ra . — 26 mV 
(D112A/V116D; blue diamonds). 
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Figure 6. The D112V/V116E mutant is a functional pro- 
ton-selective channel. (A) Families of currents generated 
by D112V/V116E in a HEK-293 cell in 10-mV increments 
at pH 0 5.5 and pli 5.5 (A), and pH 0 5.0 and pH ; 5.5 (B). 
As indicated, Vhoid was — 60 mV, and the membrane was re- 
turned to Vhoid after pulses. (C) V r „ measured at pH„ 5.0, 
5.5, 6.5, 7.0, and 7.5, and at pH; 5.5 or 6.5, indicates that 
D112V/V116E is proton selective. Inset shows tail currents 
at 10 (uppermost) , 5, and 0 mV at pH Q 5.5 and pHi 5.5. Pro- 
ton selectivity is imparted by either Asp or Glu at positions 
112 or 116, indicating that side-chain length is not critical 
for this function. 



in a cell with D112V/V116S in symmetrical pH 5.5 
TMACH 3 S0 3 solutions (Fig. 7 A), with V rev near 0 mV 
(Fig. 7 B) . Replacing bath CH 3 S03 with CI increased 
the outward current (Fig. 7 C) and shifted V rev strongly 
negatively (Fig. 7 D): on average, by —37.0 ± 2.4 mV 
(mean ± SEM; n = 7) for Asn and -35.8 ± 1 .8 mV (n = 4) 
for Ser. These values are similar to those obtained previ- 
ously for single mutants D112N and D112S upon CI 
addition: —33.1 and —40.8 mV, respectively (Musset 
et al., 2011). Thus, the introduction of Glu, Ser, or Asn 
at position 116 conferred permeability onto the non- 
conducting D112V mutant, in each case recapitulating 



the selectivity of the corresponding single mutations at 
position 112. Viewed in terms of the ability of a single- 
amino acid substituent to produce anion or proton se- 
lectivity, positions 112 and 116 were identical. 

Zn 2+ sensitivity of Arg— >His mutants shows S4 position 
in open channels 

It seemed possible that shifting the selectivity filter out- 
ward might alter the position of S4 in open channels. As 
a frame of reference, R21 1H currents were found previ- 
ously to be sensitive to internal but not external Zn 2+ 
during depolarizing pulses, suggesting that Arg 211 remains 



CH 3 SQ 3 




D112V V116S 



Figure 7. The D112V/V116S mutant is 
permeable to CP. Replacing CH 3 S0 3 ~ 
with CP increased outward current in fam- 
ilies (A and C) and shifted the tail current 
reversal potential negatively (B and D). 
All measurements were in a COS-7 cell at 
symmetrical pH 5.5 with TMA + CH 3 S0 3 ~ 
solutions or external TMA + CP solutions, 
as indicated. Pulses were applied in 10-mV 
increments from —40 to 80 mV for fami- 
lies and in 10-mV increments for tail cur- 
rents as indicated. The cell was held at 
— 30 mV and returned to Vhoid after pulses. 
The same calibration bars apply to families 
(A and C) and tail currents (B and D) . 
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internally accessible in the open state (Kulleperuma 
et al., 2013). The third Arg in S4 of the Ciona H v l was 
reported to be internally accessible in closed but not 
open channels (Gonzalez et al., 2010), perhaps reflecting 
the large size of the MTSET probe. Internal accessibil- 
ity of position 211 in the filter-shifted R211H/D112V/ 
VI 16D mutant was evaluated by testing the Zn 2+ sensitiv- 
ity of inside-out patches. This construct was proton se- 
lective (Table SI). In nine patches with R211H/D112V/ 
VI 16D, 10 pM Zn 2+ was introduced into the bath at pH ; 
7.0, and test pulses were applied. The H + current was 
inhibited by 90.0 ± 3.6% (mean + SEM), significantly 
(P < 0.0001 by unpaired t test) more than 29.7 ± 6.8% 
in six control patches (D112V/V116D). These results 
show that the 211 position is accessible to internally ap- 
plied Zn 2+ , but they do not distinguish whether the site 
is accessible in closed or open channels. 

Accessibility of His 211 in R211H/D112V/V116D chan- 
nels in the open state was evaluated by adding Zn 2+ or 
EGTA during pulses. In Fig. 8 A, shortiy after the start of 
a pulse, 10 pM Zn 2+ was introduced (red record), pro- 
ducing slowly progressing block. The subsequent pulse 
(Fig. 8 A, blue) illustrates the full extent of inhibition. 
Given that the probability of being open is high dur- 
ing large pulses — 95% at pH ; 5.5 and 75% at pH ; 6.5 
(Cherny et al., 2003) — and that gating is slow, cycling 
of channels through closed states during the pulse 
seems unlikely. Because Zn 2+ must diffuse through the 
unstirred volume at the tip of the pipette to reach the 
membrane, onset of block was slow and likely dependent 
on pipette and patch geometry, but was observed un- 
equivocally in four patches during long pulses. The ad- 
dition of EGTA during a pulse to remove Zn 2+ resulted 
in rapid reversal of block (Fig. 8 B) in six patches. The 
application of 10 pM Zn 2+ to inside-out patches from 
cells expressing the control construct D112V/V116D 
produced comparatively minor effects on the current 
(Fig. 8 C), reminiscent of its effects on native H v l 



(Cherny and DeCoursey, 1999). Thus, His at position 
211 in R211H/D112V/V116D was accessible to the in- 
ternal solution even in the open state, indicating similar 
accessibility of the S4 Arg residues in open channels 
with the selectivity filter at either position 112 or 116. 

Introducing Asp into S2 or S3 did not support 
proton conduction 

Given that Asp can produce proton selectivity at two lo- 
cations on SI , we moved Asp to several positions in the 
S2 or S3 transmembrane segments, always in the non- 
conducting Dl 12V background. We chose five locations 
in the outer vestibule (S143D, I146D, L147D, V178D, 
and S181D) that face the pore in our model and are 
located roughly between the levels of positions 1 12-116. 
Four failed to produce distinct current. D112V/S143D 
generated small Zn 2+ -insensitive currents that reversed 
at - 17 ± 6.6 mV (mean ± SEM; n= 5) at pH 0 7.0 and pH; 
5.5, well positive to En, which is —87 mV. In five cells, 
V [ev shifted by -24.3 ±1.6 mV when CH 3 S0 3 was re- 
placed with CT at pH Q 5.5 and pH; 5.5, demonstrating 
CI permeation. That Asp at position 143 (in the S2 seg- 
ment) overcame the nonconduction of Dl 12V is consis- 
tent with position 143 facing the pore, as predicted in 
our model, but in a location incompatible with its pro- 
ducing H + selectivity. 

MD simulations reveal significant differences 

in the electrostatic properties of mutant channels 

To assess the structural impact of the filter shift, MD 
simulations were performed on the WT protein, single- 
point mutants D112V and D112S, and double mutants 
D112V/V116D (VD) and D112V/V112S (VS) in a hy- 
drated lipid bilayer bathed in 500 mM NaCl, based on 
our homology model (Kulleperuma et al., 2013). The 
overall structure of the channel was preserved in the 
mutants. In particular, the average root-mean-square 
deviation between WT and mutants ranged from 2.0 A 
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Figure 8. The third arginine position in the S4 segment is accessible to the internal solution in the open state. (A) The black trace 
shows proton current in Dl 12V/V1 16D/R21 1H in an inside-out patch from a HEK-293 cell. 10 uM Zn 2+ was introduced into the internal 
solution during the red trace, and blue is the subsequent pulse in the continued presence of Zn 2+ . (B) The blue record is from D112V/ 
V116D/R211H in the presence of 10 pM Zn 2+ , and the red record shows recovery from block on the addition of EGTA shortly after 
the start of the pulse. Black represents after washout. (C) The control mutant, D112V/V116D, in a COS-7 cell exhibits weak IC Zn 2+ 
sensitivity: black, before; blue, in 10 pM Zn 2+ ; gray, after washout. All pulses are to +50 mV at pH[ 7.0 and pH„ 7.0 in inside-out patches 
of membrane. 
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(D112S) to 2.2 A (D112V). In all five systems, the aver- 
age axial position of helices S1-S4 relative to the rest of 
the bundle was 0.73 ± 0.06 A, 0.34 ± 0.05 A - 1 .80 ± 0. 1 6 A, 
and 0.53 ± 0.06 A (SEM), respectively, indicating that 
the registry of the four helices was conserved in all the 
simulations. Combining data from the five systems leads 
to axial distributions that fit a single Gaussian distribu- 
tion for each of the four helices, S1-S4, with standard 
deviations cr = 0.4, 0.5, 0.7, and 0.4 A, respectively, 
emphasizing the small amplitude of axial fluctuations. 
Furthermore, the average hydration profiles of the pore 
in WT and various mutants did not differ significantly 
from one another (Fig. 9). The pore is characterized by 
an hourglass shape with an ^13-A-long narrow region 
comprising an IC bottleneck at Phe 150 and an extra- 
cellular (EC) botdeneck at Asp 112 . It is noteworthy that 
average pore hydration was quite similar for proton- 
selective (WT, VD), anion-permeable (VS, S = D112S), 
and nonconducting mutants (VAL = D112V), indicating 
that the mean hydration profile is not a good predictor 
of selectivity. 



Despite these overall similarities (Fig. 9), differences 
were evident in the local structure of the pore near the 
EC bottleneck (Fig. 10). In the WT (Fig. 10 A), the EC 
constriction usually consists of a salt bridge between 
Arg 208 and Asp 112 (or occasionally Asp 185 ), whereas the 
other charged residues in the EC vestibule, namely 
Asp 185 and Arg 205 , usually form a spatially distinct ion 
pair. Consistent with our previous simulation study 
(Kulleperuma et al., 2013), Asp 11 2 -Arg 208 in WT is pres- 
ent as an ion pair most of the time but is occasionally 
disrupted by water molecules, resulting in the transient 
appearance of a water chain. In all four mutants consid- 
ered, however, the absence of a charged side chain at 
position 112 led to the reorganization of the ionic net- 
work in the EC vestibule. In the proton-selective VD mu- 
tant (Fig. 10 B) , various arrangements of ionic networks 
involving between two and four charged side chains 
from helices SI (Asp 116 ), S3 (Asp 185 ), and S4 (Arg 205 and 
Arg 208 ) were observed. Fig. 11 illustrates the most fre- 
quent configurations adopted by Asp 116 . Most of the 
time (^82%), it interacts with Arg 205 (Fig. 11 A) or Arg 208 
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Figure 9. Pore hydration is similar in WT and 
several mutant channels despite very different 
selectivity. Average water density within a 0.7-nm 
radius of the mean axis of the pore is plotted, 
normalized to the bulk water density for 5,000 
snapshots from each replica of different sys- 
tems. The membrane boundaries are indicated 
by dashed lines, with the external surface to 
the right. The nadir is near Phe 150 in all cases. 
Average axial water density for: (A) WT (proton 
selective) and D112V (VAL, nonconducting); 
(B) WT and D112V/V116S (VS) andD112S (S),two 
anion-permeable channels; (C) WT and D112V/ 
V116D (VD), of which both are proton selective. 
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(Fig. 11 D) or both (Fig. 11 C). Similarly, in WThH v l, the EC funnel does not depend on whether or not 
Asp 112 was engaged in salt linkage 90% of the time but Arg 208 is engaged in ion pairing (Fig. 10 B). Further- 
almost exclusively with Arg 208 . However, in WT hH v l, a more, our models predict the presence of water path- 
continuous water chain was present only when the ways in the anionic mutants D112S and VS (Fig. 10, 
Asp 112 -Arg 208 salt bridge was broken, which occurred C and D), and even in the nonconducting D112V mu- 
~10% of the time (Kulleperuma et al., 2013). In con- tant (Fig. 10 E). In all systems considered, the channel 
trast, water pathways were observed in all configura- contains a narrow bottleneck M). 5 nm in length ( — 0.68 
tions of the VD mutant, reflecting the greater width of < z < —0.18 nm) between the tips of the EC and IC fun- 
the pore at this level compared with position 112. That nels. Although this botdeneck is lined with nonpolar 
an aqueous pathway is not predictive of proton selectiv- residues including Val 109 , Phe 150 , Val 178 , and Val 209 , it con- 
ity is not surprising given the example of aquaporin tains a single file of two water molecules at least 75% of 
channels that conduct water at high rate but are imper- the time, consistent with a putative ion pathway between 
meable to protons (Zeidel et al., 1994). In contrast with the EC and IC funnels. 

WT, pore hydration is not significantly modulated by As a first step toward characterizing the energetics of 

the configurations of Asp 116 in the VD mutant. ion permeation in the channel, we computed the static- 

In particular, Arg 208 was observed to form ion pairs field energy for a virtual positive point charge along water 

with Asp 116 , Asp 185 , or neither, or both. In contrast to WT, pathways spanning the length of the channel. Fig. 12 A 

the local hydration of the constriction at the bottom of shows that D112V, D112S, and VS mutants contain an 




Figure 10. The EC salt-link network realigns in mutants. Representative snapshots of ionic networks in the WT and mutants, with the 
external end up. (A) WT protein in contact (left) and water-mediated (right) states of the D112-R208 ion pair. (B) VD mutant, with 
R208 participating in (left) and free from (right) the EC salt-link network. (C) VS mutant. (D) D112S mutant. (E) D 1 1 2V mutant. Acidic 
and basic side chains in the external and internal funnel, together with side chains of residues 112, 116, F150, and R211 are shown in 
licorice representation together with ribbon traces of the four a-helical transmembrane regions: red, SI; yellow, S2; green, S3; blue, S4. 
Salt bridges are shown as orange lines. The volumetric surface of water within the pore is shown with a water radius of 0.14 nm. 
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electrostatic barrier opposing the movement of a cation 
through the EC side of the bottleneck, which is lacking 
in the WT. In contrast, the energetic properties of the 
VD double mutant depend on the local arrangement of 
ionic residues in the EC funnel. Specifically, a static- 
field barrier to cation movement is present at the EC 
botdeneck when Arg 208 does not take part in any salt 
link (Fig. 12 A, "VD unpaired"), but this barrier is re- 
duced by at least one half when the guanidinium group 
of Arg 208 is paired with the carboxylate group of either 
Asp 116 or Asp 185 , or both (Fig. 12 B, "SL," salt linkage). 
In the latter cases, the static-field profile becomes com- 
parable to that of WT. Intriguingly, the electrostatic pro- 
file seems to be relatively insensitive to the nature of the 
ionic pairing of Asp 116 (illustrated in Fig. 11). 

Collectively, the above results suggest that the lack of 
proton selectivity in D112V, D112S, and VS mutants is 
caused at least in part by the distribution of charged 
groups in the lumen, where the presence of an excess 
positive charge near the EC constriction would tend to 
favor anions over cations. This finding is consistent with 
the anionic selectivity of D112S and VS mutants. In- 
versely, the neutrality of the EC botdeneck region in 
WT and VD systems results in approximate cancellation 




Figure 1 1 . Several configurations of Asp 116 in the Dl 12V/V1 16D 
mutant. Two-dimensional histogram of the distances from the 
center of charge of D116 to that of R205 (di) and R208 (<k), re- 
spectively, with increasing probability of ion pairing from purple 
to green. Snapshots surrounding the graph illustrate each type 
of interaction circled in the graph: (A) linked D116-R205 pair, 
37% of the time; (B) open (no salt bridge), 18%; (C) D116 linked 
to R205 and R208 simultaneously, 33%; (D) D116-R208 salt 
bridge, 12%. 



of the static field, which, although it does not explain 
proton selectivity, is consistent with the fact that these 
two channels are permeable to a cation. 

DISCUSSION 

Nonconducting mutants 

When the critical Asp 112 residue was moved along the 
SI segment from position 108 through 118 (withD112A) , 
8 of 10 mutants did not conduct, proton current was 
seen at 1 16, and anion current was seen at 109. Positions 
116, 112, and 109 all face the "pore" in the predicted 
open-state structure of hH v l (Kulleperuma et al., 2013). 
Asp at pore-facing position 143 in the S2 segment also 
exhibited anion current. Evidently, the Asp carboxyl 
group must face the pore to enable conductance of any 
kind, and with Asp at a nonpore-facing location, the SI 
segment is not free to rotate enough for the carboxyl 
group to reach the pore. When Asp faces away from the 
pore, its pK^ likely increases substantially, making it per- 
manently neutral. When Lys was introduced at a series of 
locations in the acetylcholine receptor channel, its pK^ 
was decreased, often drastically, when it did not face the 
pore directly (Cymes et al., 2005) . When ionizable amino 
acids are inserted inside proteins by mutation, their pK^ 
generally shifts in the direction that promotes neutrality 
(Isom et al., 2008). In contrast, in native proteins, ioniz- 
able residues have evolved to establish interactions with 
their neighbors that favor ionization (Kim et al., 2005; 
Gunner etal., 2011). 

The rules of the game 

The selectivity of various mutants summarized in Table 1 
reveals "the rules of the game" for hH v l. In the mutants 
studied, the presence of Asp or Glu at position 1 1 2 or 1 1 6 
was necessary and sufficient to produce H + selectivity. 
The identity of the amino acid at position 109 had no 
effect; selectivity was determined entirely by positions 
112 and 116. When Val was present at both critical posi- 
tions, 112 and 116, the channel failed to conduct (an 
effect that was not overcome by introducing Asp at 109) . 
That Val seems to disfavor permeation may reflect its 
relative hydrophobicity (Kyte and Doolittle, 1982; Hessa 
etal., 2005). The VxxDxxxV motif (or a closely conserved 
version in which Leu replaces Val 109 and/or Val 116 ) ap- 
pears in the SI segment of all eight species with electro- 
physiologically confirmed H v l (DeCoursey, 2013). In 
hHvl, when positions 112 and 116 were occupied by Val 
and a small neutral amino acid (Ala or Ser, for example) 
in either order, the result was anion permeation. 

Microenvironment of positions 112 and 116 
With respect to selectivity, voltage dependence, kinet- 
ics, and ApH dependence of gating, the filter-shifted 
D112V/V116D channels were quite similar to WT. In 
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Figure 12. Effect of the charge 
distribution of the channel on 
the energetics of ion transloca- 
tion in WT and mutant chan- 
nels. Orientation is as in Fig. 9, 
with the external solution to 
the right. Static-field energy for 
the transfer of a positive point 
charge in (A) WT, D112V, VD 
when R208 is unbound to the EC 
salt-link network, and two an- 
ion-permeable channels, VS and 
D112S. (B) WT and VD when 



-3-2-1 0 1 2 3 R208 is forming salt link(s) (SL) 

z (nm) with VI 16D or D 185, or both. 



addition, the effects of point mutations at position 116 
closely resembled those at 112. Glu replacing Asp at 
either position preserves proton specificity. Channels 
with neutral residues like Ser or Asn at either position 
(with Val occupying the other) conduct anions. These 
phenomenological parallels indicate that in terms of 
the molecular details required to establish proton (or 
anion) selectivity, the two positions are virtually indis- 
tinguishable. However, the identities of individual 
neighbors of the critical aspartate differ substantially 



at the two positions. Side chains with atoms within 6 A 
of Asp 112 (WT) or Asp 116 (in the Dl 12V background) at 
least 50% of the time were identified from time-averag- 
ing of MD simulations. Intriguingly, in WT, the nearest 
neighbors to Asp 112 , excluding those on the SI segment 
(Ala 113 , Leu 111 , Leu 115 , Val 116 , and Val 109 ), include resi- 
dues on both S2 (He 146 , Ser 143 , and Phe 150 ) and S4 (Arg 208 ) , 
whereas in the VD mutant, only S4 residues (Arg 205 and 
Arg 208 ; along with SI residues Leu 117 , Glu 119 , Leu 115 , and 
Leu 120 ) are within 6 A of Asp 116 . In the D112V/V116D 

TABLE 1 



Effects on 


selectivity of amino 


acids at positions 109, 112, and 116 in 


the SI transmembrane segment ofhHyl 


Construct 


109 


112 


116 


Selectivity 


Reference 


WT hH v l 


Val 


Asp 


Val 


H + 


Many 


A113D 


Val 


Asp 


Val 


H h 


This 


D112E 


Val 


Glu 


Val 


H h 


Musset et al., 2011 


D112A/V116D 


Val 


Ala 


Asp 


H h 


This 


D112V/V116D 


Val 


Val 


Asp 


H h 


This 


D112V/V116E 


Val 


Val 


Glu 


Ff 


This 


D112V 


Val 


Val 


Val 


0 


Musset et al., 2011 


D112V/V109D 


Asp 


Val 


Val 


0 


This 


D112A, N, S, H, K, F 


Val 


Ala 


Val 


er 


Musset etal., 2011 


D112A/V109D 


Asp 


Ala 


Val 


er 


This 


D112V/V116S 


Val 


Val 


Ser 


er 


This 


D112V/V116N 


Val 


Val 


Asn 


cr 


This 



Selectivity to H + means V^ ev was close to £ H at various ApH; 0 means no credible currents; CI means l^ev shifted negatively when CI replaced CH 3 S0 3 
the external solution. Acidic residues are bold; neutral residues other than Val are italicized. 
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mutant, the selectivity filter is roughly twice as far from 
Phe 150 (the external delimiter of the charge transfer 
center; Tao et al., 2010) as Asp 112 is in the WT channel. 
This result means that immediate proximity of the car- 
boxyl group to the most hydrophobic point in the chan- 
nel, which in all constructs occurs at Phe 150 (Figs. 9 and 10), 
is not essential to establishing any of the main features 
of the channel, including proton specificity, gating ki- 
netics, and ApH-dependent gating. 

Why does 109 not work? 

The lack of influence of position 109 (Table 1) suggests 
that selectivity is established in the external vestibule, 
outside the charge transfer center boundary at Phe 150 . 
All current models of hH v l place Asp 112 in the external 
vestibule (Musset et al., 2010; Ramsey et al., 2010; Wood 
et al., 2012; Kulleperuma et al., 2013). The MD simula- 
tions provide additional insight. In WT hH v l, Asp 112 is 
frequendy paired with Arg 208 . During the 10% of the 
time that these residues were unpaired, a continuous 
water wire was present. With Asp moved to 116, there 
was an electrostatic barrier to cation permeation when 
Arg 208 was unpaired, but not when it was paired with ei- 
ther Asp 185 or Asp 116 . In support of these findings, an 
unpaired Arg residue near the EC mouth of aquaporin 
channels is essential to block the translocation of pro- 
tons (Beitz et al., 2006; Wu et al., 2009) . 

Limited plasticity of the channel with respect to axial helix 
movement does not support a change of S1 registry 
in the mutants 

The interpretation of measurable effects of point muta- 
tions assumes conservation of the global protein con- 
formation. Accordingly, we have modeled the structure 
of the mutants by substituting individual side chains in 
our model of the WT. An alternative interpretation for 
the native-like phenotype of the VD double mutant is 
that the SI helix is shifted by one turn toward the inte- 
rior of the cell, so that Dl 12V and V116D take the place 
of VI 09 and D112, respectively. However, the analysis 
of collective fluctuations from the simulations shows 
that not only were the overall structure and hydration of 
the channel conserved, but there was also no significant 
axial displacement of individual helices in response to 
any of the mutations, and the plasticity of the channel 
with respect to axial helix movement remained limited 
in all the systems studied (Figs. SI and S2) . Although 
the possibility of a change in SI registry with respect to 
the rest of the bundle cannot be ruled out, these results 
suggest that such a change is unlikely in any of the mu- 
tants considered in the present study. 

Collectively, our current findings suggest that charge 
neutrality resulting from the presence of Asp (or Glu) 
vis-a-vis Arg 208 at the external constriction is required 
for the charge selectivity of the channel; in particular, 
the presence of an excess positive charge at the narrow 



end of the external funnel in neutral Asp 112 mutants 
leads to a barrier opposing cation movement but com- 
patible with anion selectivity. In addition, the presence 
of an ionic network involving Arg 208 results in dynamic 
fluctuations of pore hydration and/ or electrostatic pro- 
perties that may contribute to the mechanism of proton 
selectivity in the WT and VD channels. 

Parallels in other molecules 

To our knowledge, hH v l is the first example of a selec- 
tive ion channel whose selectivity filter can be moved by 
a pair of point mutations. The acetylcholine receptor 
channel remains nonselective among cations when its 
ring of Glu residues is shifted by one turn of the a helix 
(Cymes and Grosman, 2012). Several other proton-con- 
ducting pathways permit shifts of critical amino acids, 
although in these molecules the portable function is 
primarily rapid proton flux rather than proton selectivity 
per se. For example, in FiF 0 -type ATP synthase (Esche- 
richia coli) , proton translocation is preserved when Asp 61 
is shifted to position 24 on another helix (Miller et al., 
1990). Asp 61 can be replaced by Glu 61 , but with dimin- 
ished proton transport, suggesting that precise location 
of the carboxyl is critical (Fillingame, 1990). Portability 
of Asp shows that the precise structure of apolar neigh- 
bors of the carboxyl group is not critical (Fillingame, 
1990). Another essential residue in ATP synthase is 
Arg 210 , which is thought to lower the pi£ a of Asp 61 tran- 
siently to ensure proton release (Fillingame et al., 2003; 
von Ballmoos et al., 2009) . 

As in ATP synthase, an Asp 213 critical for proton trans- 
location in nicotinamide nucleotide transhydrogenase 
can be replaced by Glu, but activity is decreased to 18% 
(Yamaguchi et al., 2002) . In enzyme studies, activity that 
is rate-limited by proton translocation is assessed. The 
correlate in hH v l is single-channel conductance, which 
was not examined here. Our criterion is perfect proton 
selectivity, which is preserved when Glu replaces Asp 
at either position 112 or 116 in hH v l. We cannot say 
whether Glu is equally efficient. 

When Asp 135 in the proton entry channel in cyto- 
chrome bos, ubiquinol oxidase of E. coli is neutralized, its 
function can be restored by shifting Asp to position 139 
or 142 (Garcia-Horsman etal., 1995). 

Neutralizing Asp 132 , the namesake of the D channel 
in cytochrome c oxidase from Rhodobacter sphaeroides, by 
mutations D132N or D132A nearly abolishes proton 
uptake (Fetter et al., 1995), which is restored by reposi- 
tioning the Asp at N139D (Varanasi and Hosier, 2011, 
2012). Intriguingly, proton uptake is also restored by 
removing subunit III (Adelroth and Hosier, 2006) and 
also in the D132N/N139T double mutant, showing that 
in the D channel, rapid proton uptake can be accom- 
plished without an acidic group, although enzyme turn- 
over remains impaired (Johansson et al., 2013). Asn 139 is 
thought to serve a special "gating" function in cytochrome 
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c oxidase that may normally limit WT H + flux (Henry 
et al., 2009); Thr 139 appears to optimize aqueous connec- 
tivity within the pore (Johansson et al., 2013) . 

Finally, His 64 shuttles protons from the catalytic center 
of human carbonic anhydrase II, and function is pre- 
served with His shifted to His 67 (H64A/N67H), but not 
His 62 (H64A/N62H) , despite crystal structures indicat- 
ing that the side chains of both His 62 and His 67 extend 
into the active-site cavity at distances from the catalytic 
zinc similar to His 64 (Fisher et al., 2005). 

These examples show that an amino acid side chain 
must be positioned correcdy to maintain a high rate of 
proton transfer, but that in some cases a reasonable rate 
of proton transfer can be retained upon moving the 
side chain, especially if it is not moved too far. Given 
that proton transfer via titratable amino acid side chains 
is a way to achieve proton selectivity, proton transfer and 
proton selectivity may be the same process, and thus in 
some cases, a side chain essential for H + selectivity may 
be moved without losing selectivity. 

In summary, shifting Asp along the SI segment iden- 
tified three locations in hH v l that line the pore and 
permit conduction: 109, 112, and 116. Asp produced 
proton specificity only at positions 112 and 116. When 
introduced at nonpore-facing positions, Asp abolished 
function. Glu at either position preserved selectivity, in- 
dicating leeway in side-chain length. We conclude that the 
minimal requirements for proton specificity of hH v l 
include Asp or Glu, which must face the pore directly, 
and evidendy must be located in the external vestibule, 
above the charge transfer center, Phe 150 . The portabil- 
ity of the selectivity filter indicates latitude in the requisite 
local environment. This observation seems consistent with 
the suggestion that ionizable groups that enable proton 
transport could have evolved by random mutation without 
the need to simultaneously develop a specialized micro- 
environment for charge stabilization (Isom et al., 2008) . 
On the other hand, the inability of Asp to produce H + - 
selective conductance at most positions tested argues that 
additional factors are involved. Modeling indicates fre- 
quent salt-bridge formation between Asp and several part- 
ners, which may contribute to selectivity, proton transfer, 
and conformational stability. Ionizable residues favor the 
incorporation of water molecules (Pless et al., 2011) that 
may also facilitate proton translocation. However, neither 
the presence of a continuous water pathway nor the mean 
hydration profile was found to have any clear relationship 
with the selectivity of the channel. Detailed computational 
studies of H + and other ions in the permeation pathway 
will be required to examine the molecular basis of ion 
movement and achieve a full understanding of the mecha- 
nism of proton selectivity in hHvl. 
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